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Abstract 



The consequences of adding random perturbations (anarchy) to a baseline hier- 
archical model of quark masses and mixings are explored. Even small perturbations 
of the order of 5% of the smallest non-zero element can already give deviations sig- 
nificantly affecting parameters of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, 
so any process generating the anarchy should in general be limited to this order of 
magnitude. The regularities of quark masses and mixings thus appear to be far from 
a generic feature of randomness in the mass matrices, and more likely indicate an 
underlying order. 
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1 Introduction 



The origin of fermion masses and mixings is one of most important issues in parti- 
cle physics. Unfortunately, these parameters are inputs in the well-tested Standard 
Model. All one can do is to measure them as accurately as possible and hope that 
in the future a more fundamental theory will actually be able to predict their values. 
Many attempts have been made in this direction, with grand unified theories (whether 
super symmetric or not) being favorite candidates |IJ]. However, such a fundamental 
theory could well be quite complicated, with many new fields and couplings. Conse- 
quently, at low energies, the observed fermion masses and mixings could be the result 
of a large number of contributions. If that were the case, one might expect that at low 
energy scales the mass matrices of fermions would have a random nature. This idea 
was first suggested by Froggatt and Nielsen 0, who performed a statistical study of 
fermion masses without success. 

More recently, the idea of flavor anarchy was introduced in order to explain new 
data on neutrino masses and mixings || , and several analyses based on this idea have 
been performed [Q. It was suggested that a similar model could also explain masses 
and mixings in the quark and charged lepton sectors. The purpose of this note is to 
study the robustness of such a model for random quark matrices in the quark sector. 
We take a baseline model which approximately reproduces the observed masses and 
mixings and study its sensitivity to random perturbations of the parameters. In this 
way we determine to what extent the observed regularities can survive effects which 
may be purely coincidental or generic. We find that quark masses and mixings appear 
to be far from a generic feature of randomness in the mass matrices, and more likely 
point to an underlying order. 

2 A simple ansatz 

The quark and charged lepton sectors are fundamentally different from the neutrino 
sector because of the existence of a large mass hierarchy between the families and the 
resulting small mixing. We shall construct a "baseline" description of quark masses 
and mixings which incorporates several approximate regularities. For this purpose 
we begin with quark masses evolved via the renormalization group to a common high 
mass scale \i = Mz- At this scale, the masses have been found to lie in the range || 
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These masses thus are compatible with the hierarchy 
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which we shall incorporate into our ansatz for quark mass matrices. (This is certainly 
not a property of the charged leptons, for which B 
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is a better approximation.) For illustrative purposes we will adopt e 
and edown = £d = 0.21, which approximately reproduces the observed hierarchies. 
We also seek a set of mass matrices reproducing the regularities 
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The first of these was noted some time ago J?], ||. We further wish to reproduce 
the hierarchy of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements noted by 
Wolfenstein [IIJ, in which 
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(A ~ 0.22), and with |K&| < |Kd| as favored by fits to data []TT] . 

These regularities can be incorporated into a simple quark mass ansatz in a basis 
which we call hierarchical: 
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where m3 denotes the mass eigenvalue of the third-family quark (t or 6) . Hierarchical 
descriptions of this type were first introduced by Froggatt and Nielsen [0. The 
present ansatz is closely related to one described by Fritzsch and Xing [ 131, IHI- We 
shall not be concerned with relative coefficients of order 1 in different terms; for 
example, models in which the off-diagonal e 2 terms are multiplied by v^2 may fit 
| Kb | somewhat better [|l^, 15]. The e 3 terms in our model are separate parameters 
in Ref. [Q. We have assumed the phase to be present only in the e 3 terms; there is 
little sensitivity to phases in the off-diagonal e 2 terms |T^]. Our purpose is primarily 



2 



to construct an easily-manipulated "cartoon" version of the mass matrices, so as 
to study the robustness of their predictions for masses and mixings under random 
perturbations. 

The eigenvalues of the matrix (^) to order e 4 are given by 

Ai = -m 3 e 4 (8) 
A 2 = m 3 e 2 (9) 
A 3 =m 3 (l + e 4 ) (10) 

and are independent of the phase <fi. Therefore, in this ansatz the quark masses 
naturally obey the hierarchy (|3|). 

The matrix ([?]) can be made real by a unitary transformation M.^ = PAdnP^, 
where P is given by 

(1 0\ 
p = . (11) 

\0 e**J 

The real symmetric matrix then can be diagonalized by an orthogonal trans- 
formation T M^O where the matrix O is given to order e 4 (its columns are u,c,t) 
by: 

O = I + f 1-4-T ^ + f I • (12) 
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2 

We will assume this same mass matrix ansatz for both the up and the down 
quark sectors. In this case, the Cabibbo-Kobayashi-Maskawa (CKM) matrix is readily 
obtained by: 

VcKM = 0^ p P U pPj own Odown • (13) 

We first give approximate expressions for the CKM matrix elements: 

(14) 
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where A = <p u — <p<±. These can be brought into a form closer to the standard phase 
convention (see, e.g., fllUj) by multiplying the c and t rows by e^ x_A ) and the s and b 
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columns by e lx , where x — Arg(ed — e iA e u ) is chosen so as to make V us and V C( i real. 
Then we find, to leading order in small terms, 



Kd = l-(e 2 u + e 2 d )/2 + e lA e u e d (23) 

Ks = |e d -e lA e u | (24) 

Kb = e^-^-e 2 .) (25) 

Kd = -|e d -e iA e u | (26) 

V CB = l-(e 2 +e 2 )/2 + e - iA e u e d (27) 

Kb = e 2 d -tl (28) 

Kd = e d ^(el-el) (29) 

Ks = e*-e 2 (30) 

V th = 1. (31) 



The angles in the unitarity triangle can be expressed very simply in terms of these 
quantities. We find 



a(= 02 ) = A , (3(= 0i) = -x = tan 1 



sin A 



v £d/ e u — cos A y 

7 (=0 3 ) = 7r-a-/3 . (32) 
These expressions also hold in more general versions of the present model [O, [L4j]. 



Note that \V u b\ and |T4d| are specified entirely in terms of the e Uj d, with 

\V ch \ = \V U \ = e 2 - e 2 = 0.0392 , \V ub /V td \ = e u /e d = 1/3 (33) 

for our choice of parameters. The shape of the unitarity triangle is determined entirely 
by the magnitude of V us , which changes as A varies. It has been noticed previously 
that a value of A close to 90° gives a good fit to \V US \ [|, ||. For A = 90° we find 
\V US \ = \Vcd\ = (e d + e 2 ) 1//2 = 0.221. The Wolfenstein parameters p and 77 defined by 



^p 2 + v 2 =fe} =e " ( 34 ) 

A^Vc 2 - p) 2 + V 2 = 14 = (35) 



where A = |V^g| and c = 1 — A 2 /2, are given by: 
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With our choice of parameters, p = 0.12 and rj = 0.29. 



V = \I^-P 2 - (37) 
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We now diagonalize the u and d quark mass matrices exactly (without the above 
approximations), and calculate the resulting CKM matrix. As an illustration, for the 
same parameters as above, we find 



|Kd| = \V a \ = 0.975 (38) 

|Ks| = |Kd| = 0.222 (39) 

I Kb | = |Ks|= 0.0410 (40) 

V td \ = 0.0080; |Kb| = 0.0029; | Kb/Kb I = 0.071 (41) 

a = 98° ; (3 = 18° [sin(2/3) = 0.58] ; 7 = 64° (42) 

p = 0.15 ; 7] = 0.29 . (43) 



These predictions are reasonable and we will consider this simple model as our 
baseline model. We now perform random perturbations around it (anarchy) to study 
the robustness of these predictions. 



3 Anarchy 

The underlying theory of fermion masses is unlikely to guarantee that a particular 
mass matrix element is absolutely zero or has a fixed value given by high energy 
symmetries of the theory. One could expect radiative effects or small symmetry 
breaking parameters to contribute to a given mass matrix texture. Since we don't 
have a complete theory to compute them, we will explore the consequences of adding 
randomly small perturbations to our baseline model. These perturbations are what 
we call anarchy. We will study just how much anarchy can our baseline model tolerate. 

We will add to each element of the up and down mass matrices a random variable 
£ with real part given by a random gaussian of zero average and a phase given by a 
random variable uniformly distributed in the range [0, 27r]. We ensure that the mass 
matrices remain hermitian by requiring that Qj = Q { . 

The standard deviation of the real part of these random variables will be our 
measure of anarchy. In fact, it would be natural to have the standard deviations 
proportional to the smallest non-zero elements of the mass matrix, namely a = ae 3 , 
where a is the standard deviation of the gaussian distribution. 

In Figs. 0, H], and ||] we illustrate the effect of anarchy in the p — rj plane and the 
iill _ 1112. planes for both up and down sectors for different values of a. As expected, 



mi m,2 

the dispersion in the figures increases with a. The quantities p, rj and ^ are more 
stable under perturbations than the ratio of the two lightest quark masses, 

In order to measure the degree of variation of these parameters arising from anar- 
chy, we introduce the quantity k, the ratio of the standard deviation ax to the mean 
Px for a given parameter A in a simulation: k = ax/ Px- in Table [1] we show k for 
the different parameters generated by 1000 simulations. 

We can see that k is roughly proportional to the parameter a determining the 
width of the gaussian distribution used to generate the different perturbations. The 
sensitivity in the parameters increases in the following order: (jj^j , (^f) d 5 
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Figure 1: Effect of anarchy in the p — 77 plane for different values of a. (a) a = 0.01; 
(b) a = 0.05; (c) a = 0.1. 

(^2) ' P anc ^ (^2) ' sur P r i s i n gly! ^e most sensitive parameter is (j^j ? 
since it involves the ratio of two small quantities. Values of a < 0.1 are required in 
order not to disturb the hierarchy significantly. These correspond roughly to C(e|) ~ 
0(X 5 ) entries in the mass matrices (0). 

4 Conclusions 

We have performed a simple exercise of exploring the consequences of adding ran- 
dom perturbations (anarchy) to a baseline hierarchical model of quark masses and 
mixings. We find that even small perturbations, with gaussian distribution of zero 
mean and standard deviations of the order of 5% of the smallest non-zero element 
can already give deviations of 26% in the p parameter, for instance. Therefore, we 
conclude that any physics process generating the anarchy, be it radiative corrections 
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Figure 2: Effect of anarchy in the up sector in the — — — plane for different values 
of a. (a) a = 0.01; (b) a = 0.05; (c) a = 0.1. 



Table 1: Ratios k of the standard deviations to the means for parameters generated 
by 1000 simulations. 



a = 0.01 a = 0.05 a = 0.1 

p 0.054 0.26 0.57 

r] 0.029 0.14 0.28 

(**) 0.00072 0.0034 0.0071 

2*0 0.14 0.60 0.73 
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(2*) 0.0021 0.010 0.021 

M 0.046 0.24 0.45 



7 



C\2 



0.055 



0.050 



0.045 



0.040 



0.035 



1 1 1 1 1 1 1 
1 (a) 


1 1 1 1 1 1 1 1 1 1 


L mi 




+ii5jgjflB 




— 




1 1 1 1 1 1 1 1 1 1 1 1 1 


+ + 

1 1 1 1 1 1 1 1 1 1 



0.0436 



0.0442 



m 2 /m 3 



0.10 
0.08 

02 


= (b) 


3 0.06 




^ 0.04 




0.02 




0.00 





n — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 




i i i_i i_i i i_i lj 




0.042 



0.046 



0.042 



0.046 



m 2 /ni3 



m 2 /m 3 



Figure 3: Effect of anarchy in the down sector in the — — — plane for different 
values of a. (a) a = 0.01; (b) a = 0.05; (c) a = 0.1. 
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or small symmetry breaking parameters, should be in general limited to this order 
of magnitude, unless some spurious cancellations occur. The regularities of quark 
masses and mixings thus appear to be far from a generic feature of randomness in 
the mass matrices, and probably indicative of an underlying order. 
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